THE HIGGS SECTOR AND ELECTRON ELECTRIC DIPOLE MOMENT 
IN NEXT TO MINIMAL SUPERSYMMETRY 
WITH EXPLICIT CP VIOLATION 



MUGE BOZ 
Physics Department, Hacettepe University 
Ankara, 06532, Turkey 

Abstract 

We study the explicit CP violation of the Higgs sector in the next-to-minimal 
supersymmetric model with a gauge singlet Higgs field. Our numerical pre- 
dictions show that electric dipole moment of electron lies around the present 
experimental upper limits. The mass of the lightest Higgs boson is quite sen- 
sitive to the CP violating phases in the theory. It is observed that as the 
vacuum expectation value of the singlet gets higher values, CP violation in- 
creases. 



1. Introduction 

The present bounds on the electric dipole moments (EDM) of the particles 
generate serious hierarchy problems concerning the amount of CP violation. For 
instance, the neutron EDM ^, induced by the the QCD vacuum angle (Oqcd), is 
approximately ten orders of magnitude larger than the existing bound, which is 
the source of the so-called strong CP problem. This naturalness problem has been 
solved by the Peccei-Quinn mechanism ^ which promotes Oqcd to a dynamical 
variable via the phases of the quarks or additional color triplets 2. 3, 4 

In the supersymmetric (SUSY) extensions of the standard electroweak theory 
(SM) this hierarchy problem still persists. Moreover, there appear novel sources of 
CP violation coming from the soft supersymmetry breaking mass terms. Though 
the phases of the soft terms have been shown to relax to CP-conserving points in 
the minimal model (MSSM) ^, this is not necessarily true in the non-minimal model 
(NMSSM) ^ containing a singlet. The La grangian of the MSSM consists of various 
mass parameters which are not necassarily real ^ . The phases of these parameters 
contribute to known CP violation observables such as the electric dipole moments of 
the electron, neutron, atoms, and molecules 8.940, 11, 12, 13^ Higgs system I4,i5,i6^ 
and the decays and mixings of mesons 

In addition to these CP hierarchy problems, in minimal SUSY model there is 
another hierarchy problem concerning the Higgsino Dirac mass parameter (/i) , that 
is, this mass parameter follows from the superpotential of the model and there is 
no telling of at what scale (ranging from AIw to Mpi) it is stabilized. 

The next-to-minimal supersymmetric model (NMSSM) is the most economic 
extension of the MSSM in which the /i parameter is induced by the vacuum expec- 
tation value (VEV) of an additional gauge singlet. The NMSSM not only solves the 
fj, problem by means of the VEV of the singlet, but also offers a rich phenomenology 
for colliders 
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In the NMSSM, the Higgs seetor of the MSSM (inchiding two CP-even Higgs 
bosons, a CP odd Higgs boson, and a charged Higgs pair) is extended so that there 
are three CP-even Higgs bosons, two CP-odd Higgs bosons, and a charged Higgs 
pair. Therefore, the Higgs phenomenology of the NMSSM significantly differ from 
that of MSSM As a result of the SM searches at LEP, the lower bound on 

the lightest Higgs boson mass is 115 GeV (and correspondingly tan/3 > 3.5) 
Therefore, from the searches at LEP2, the lower limit on the mass of the SM Higgs 
boson excludes the substantial part of the MSSM parameter space particularly (for 
rrit = 175 GeV) at low tan/3 (tan/3 < 3.5). However, in the low tan/? regime, 
NMSSM is in much better shape phenomenologically, since the Higgs boson masses 
are larger, and the fine-tuning is less 22,23 xhus, the phenomenological implications 
of the NMSSM can offer new opportunities at future colliders 

In the MSSM, since the tree level vacua are CP conserving, spontaneous CP 
violation (SCPV) can only occur if the radiative corrections are taken into account 
leading to a very light Higgs boson which has been discarded by the experiment 
In the NMSSM, spontaneous CP violation can occur even at the tree level, for those 
models without a discrete Z3 symmetry 25,26_ -pox models which has symmetry, 
the NMSSM can not produce SCPV at the tree level SCPV can only occur 
if the stop quark masses at the one- loop effective potential are non-degenerate 
However, at the tree level, unlike the MSSM, explicit CP violation is possible in the 
NMSSM 29.30 _ 

In the literature, the NMSSM has been studied through many researches includ- 
ing the effects of one- loop radiative corrections due to various particles, and their 
superpartners 3i>32,33,34,35^ ^j^^^ analysis of the Higgs potential with arbitrary 
number of Higgs singlets In the earlier studies, the corrections to mass matrices 
of the top and stop quarks have been calculated with either the one-loop effective 
potential or the renormalization group approach Among recent works the au- 
thors of Ref. [32] particularly addressed the effects of the explicit CP violation with 
the emphasis on the charged Higgs boson in the NMSSM, using the effective po- 
tential method, for which case the CP violating phases are induced from the stop 
and sbottom quark masses. This analysis is extended by including the contribution 
of the charginos and of the neutralinos "^^ to the one-loop effective potential, 
where the tree level CP violating phase is chosen to be equal to the one at the 
loop-level, and the complex phases of the chargino and neutralino contributions in 
the radiative corrections are taken into consideration. 

The radiatively induced CP violation effects in the NMSSM has also been ana- 
lyzed with renormalization group improvement It has been shown in Ref. [37] 
that the renormalization group analysis of radiative symmetry breaking always leads 
to tan^ values that are larger than 1. Indeed this result is a general feature of all 
phenomenologically viable theries that incorporate low energy supersymmetry and 
radiative breaking of the electroweak supersymmetry. More recent analysis on the 
Higgs mass spectrum motivated by the renormalization group analysis shows that 
to prevent the mass splitting between the light and heavy Higgs bosons from be- 
coming too large, the value of tan/3 should be kept moderate (;6 10) which we 
shall also assume in this work. 

In this work our aim is to investigate the CP violation capability of the NMSSM. 
Therefore, we limit our analysis to the effective potential with no renormalization 
group improvement. This accuracy has proven sufficient in obtaining the observable 
effects of explicit CP violation on the Higgs masses In doing this, we will keep the 
value of tan /3 moderate (< 10), as favoured by the renormalization group analysis 
of [38], which is a well- motivated parameter regime for the model under concern 
(e.g, tan/3 = 1.5 is allowed within the framework of NMSSM ^''). For the purpose of 
showing the mass spectra, we shall take tan/3 = 2. To give feeling of the sensitivity 
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of the mass on tan /3, wo shall also consider tan /3 = 10 regime. 

In the following we will study the radiativcly corrected Higgs masses and mix- 
ings, taking into account of the CP violation effects, in the parameter space allowed 
by the electron electric dipole moment (eEDM) constraint. We will base our calcu- 
lations to those of Ref. [32]. The main difference with the previous work springs 
from the fact that in our analysis the CP violating phases at the tree level and at 
the one-loop are not equal. Therefore, we analyze the mass spectra at the tree and 
one-loop levels, as well as the CP odd components of the Higgs boson by considering 
the effects of physical phases seperately. In the numerical analysis, wo particularly 
focus on the real and complex cases of the tree-level coupling A. Analyzing these 
specific cases for various parameter planes, we search for the impacts of the tree 
and one-loop level contributions on the Higgs sector which give the opportunity of 
comparing the results. In our analysis, we focus on the regions of the parameter 
space which are allowed by the eEDM constraint. 

The organization of the paper is as follows: In Section 2, we study the Higgs 
masses and the one-loop eEDM in the NMSSM. In Section 3, we discuss two special 
cases of A for two different values of Vs, and At at tan/3 = 2, and tan/3 = 10 regimes. 
In Section 4, we conclude the work. 



2. The Model 

The Higgs potential of the model at the tree level is given by : 

Vo = mij/^dp + mjjJHul'' + m||5|2 - [XAxHaH^S + ^AkS^ + h.c] 

+ '-{gf+gl){lH,f-\HX)\ (1) 

where the first, second and the third lines represent the soft SUSY breaking terms, 
F and D term contributions, respectively. Here, Ax, and Ak are the for the trilinear 
soft SUSY breaking terms. 

Assuming XA\, and kAk are real and positive, the phase between A and k* is 
given by: 

(pk\ = Arg[\k*] ^ ip^- ^Pk , (2) 

and this phase forms the source of CP violation at the tree level in Eq. (1). 

After electroweak breaking, the Higgs doublets and one Higgs singlet in (1) can 
be expanded as: 



^d \ _ f Vd + (t>i+ Wi 



and 



S = Vg + + iifs . (4) 



As usual, we calculate the Higgs masses and their mixings up to one loop accu- 
racy via 

/ d"^ V \ 

A/"^ = o o , where Xi eB= {(Ai,02,Vi,<^2,0s,Vs} • (5) 
\dXidXjJo 
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Here, V = Vq + Vi-ioop is the radiativcly corrected Higgs potential and as 
we mentioned before we take into account the top-stop and bottom-sbottom loop 
corrections. 

The stop and sbottom mass-squared eigenvalues are given by: 



tl,2 



1 



m 



Ol ■ 



1 



where the stop and sbottom mass-splittings read as: 



(6) 



A? = Ji(M2 - M^)^ + ml(^Al + X^v^,{t0y + 2A,Xv,ti3Cos^xt) ■ (7) 



Here, tj3=tanf3, t^^=cot (3. For convenience, we set the soft SUSY breaking scalar- 
quark masses as Mq = = M^, and the squark trilinear couplings as At = A^. 

The stop and sbottom mass splittings depend explicitely on the total CP viola- 
tion angle ipxt between At = At,, and A 



ipxt = Arg[XAt] 



(8) 



which forms the source of CP violation at the tree and one-loop level Higgs potential 
of the NMSSM. 

The (5 X 5) dimensional Higgs mass-squared matrix can be expressed as: 



Mij = Mij + AMij . 



(9) 



Here, comes from the tree-level potential, whereas AMij from the stop and 
sbottom contributions at the one- loop level 

The elements of the mass matrix at the tree level are are given by: 



Mi2 
Mi3 

Mm 



Ml. = -3 



M22 
M23 

M24 
M25 



[mz cos P] + [niA sin /?] , 

— [mz + mji^ — 2X V \ sin/3cosp 

, 

' V 



m\ siv? f3 cos P — 2{X'^v^) cos /? - 



(A t;^)sin/3sin(^At 



= [mz sin /3] ^ -|- [toa cos /?] ^ , 



(A v^) cos ipxt sin P 



mj^ sin/3 cos /3 — 2(A t;^) sin/3 + I -r K-^ Vg)cosipkx cos/3 



Mf 



52 



Vs 



- )(A t;Jsin(pAt cos/3 
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Mas = 




M34 = 




M35 = 






(- 


M44 = 



(A wJsinv?At 



TO^ sin /? cos /3 — 3 1 ^ 1 (''^ ) cos (p^t 



sin /3 cos /3 1 sin /3 cos fimjs^ 



(A^Ug)cOS ly^At 



+ 



M45 
where 



A2 

k 



{\vs)Ak , 



2-[(A2w2)sin2/3sin(^At] 



+ 3 



sin P cos P 
{Xvs)Ak , 



sin cos /? + 3 1 - I (A v^) cos ipxt 



{\Vs)\Ax + kvs cosy At] 
sin /3 cos /3 



(10) 



(11) 



The radiative correction terms due to the stop and sbottom corrections at the 
one-loop (AA'/jj) can be found in the work of Ref. [32]. However, as mentioned 
before, we differ from Ref. [32] in the sense that, in our analysis, the CP violating 
phase at the tree-level is not equal to the one at the one-loop. 

In our analysis, we will particularly concentrate on the lightest Higgs boson, 
whose mass can be obtained by the diagonalization of the Higgs mass-squared 
matrix by the similarity transformation: 



(12) 



where 7^.7?."^ = 1. Here, we define /15 to be the lightest of all five Higgs bosons . 

The mass eigenstates of the lightest Higgs boson (/15) can then be decomposed 
in terms of the basis elements as: 



/15 = X]^5i*i , 



(13) 



where $1, $2, '^3, ^4, $5 correspond respectively, (f>2, sin(3fi + cos(3(p2, 
ips components of the Higgs boson under consideration. 
From Eq. (13), we define the dimensionless quantity pi, 



100 X ITZsil 



1,2,3,4,5 



and 



(14) 



which is a measure of the percentage CP component of a given mass eigenstate 
Higgs boson. Therefore, in (14) for instance, ps and p^ are measures of the percent- 
age CP odd components of /15. 
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The main contributions to the one-loop eEDM come from the neutraUno, and 
chargino exchanges which can be expressed as: 



e / \ e / \ e 



By taking into account of the neutralino-selectron interaction, the neutrahno con- 
tribution to the eEDM can be written as: 

(v)'"^ = £T;{EEM.ejM,o,-,B(M^o,,j|, (16) 

which is very similar to that of MSSM except for the fact that the neutralino 
sector now extends to 5 x 5 mass matrix in the NMSSM, due to the presence of the 
additional gauge singlet. By convention, in Eq. (16) and in the formalism below, 
we use the short-hand notation, for the generic indices a and p, M^^p = and 

sw{cw) = sin 9w {cos 9 w), = cos/3, = tan6'vK- 

In Eq.( 16) 7]^^^. is the neutralino vertex which is given by the following expres- 
sion: 



Ik 



(17) 



where Af is the unitary matrix diagonalizing the neutralino matrix: Af^ M^o Af = 
diag (m^o,- ■ ■ ,?n^o j • The eigenstates (ei, 62) in Eq. (17) can be obtained via the 

unitary rotation: »S| Sg = diag (m|^ , mf^) . 

On the other hand, taking into account of the sneutrino-chargino interaction, 
the chargino contribution to eEDM reads as ^®: 

where U and V are the unitary matrices diagonalizing the chargino mass matrix: 
WMcV-'^ = diag(m^+, m^+). In Eqns. (16)-(18), B and A are the the loop 

functions 

We would like to note that we take into account one loop contributions to eEDM. 
It was pointed out in Ref. [40] that in certain regions of MSSM parameter space 
certain two-loop contributions can not be non-negligible. However, these two-loop 
contributions become sizeable only at high tan/3 (tan/3 > 30). In this work we 
have restricted our attention to moderate tan/3 values (tan/3 ^ 10) which is a 
well-motivated parameter regime for the model under consideration . Hence the 
two-loop eEDM will not provide significant contribution in our analysis. 

Eq. (15), possesses the sources of CP violation through 

{i) the gaugino masses 

(m2,mi) ^ (m2e''^^mle'^l) , (19) 
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which involves gaugino {(fi) and the SU(2) gaugino (1^2) phases in neutrahno and 
chargino mass matrices. 

{a) the complex selectron trilinear coupling 

(A,) ^ {A.e^'^') , (20) 

Therefore, the phases of the stop and sbottom trihnear couphngs {(pAt = fAh)i 
the phases of the gaugino masses {(pi, and (^2), the phases of the the Higgs potential, 
at tree and at the one- loop level {(pkx, and ipxt, respectively) form the CP violating 
sources in the full parameter space. In the following, we will perform a numerical 
study to determine the eEDM, the effects of the physical phases of the model on 
the mass and on the CP-odd components (pa and ps) of /15, in the parameter space 
allowed by the eEDM constraint. 

3. Numerical Analysis 

In this section we will consider various parameter planes to adress the issue of 
whether or not the various parameters would lead to a large amount of CP violation 
opportunities. As seen in the previous section, there are two physical phases which 
contribute to the Higgs boson mass matrix. Besides these physical phases, the free 
parameters appearing at the tree and one- loop levels are tan/3, A/., A\, k, Ak, Vg. 
At and Ai,. 

In our analysis, we particularly concentrate on h^, and analyze various parameter 
planes, in the parameter space allowed by the eEDM constraint. In doing this, we 
use the present experimental upper bound of the eEDM 41.42. 

de < 4.3 X 10-2^ e.cm , (21) 

A convinient way to observe the effects of the eEDM constraint, is via the 
dimensionless quantity: 

eEDM = , (22) 

which measures the fractional enhancement or suppression of the eEDM with respect 

to its experimental value. 

In the numerical analysis, for the purpose of definiteness, we set: {i) the trilinear 
couplings At = Ab [ii) Ax = Vg, and Ak = 100 GeV. We focus on the values of k, 
and A that are favored by the renormalization group equations in the NMSSM 
In doing this, we fix the k parameter to be near its fixed point value (i.e. we choose 
k = 0.63), and we focus on the allowed values of A, in the 0.12 < A < 0.82 interval. 

We first search whether or not the eEDM is consistent with the present experi- 
mental bounds in the above mentioned intervals provided that the gaugino masses 
are of C'(TeV). Namely, we take: {i) the gaugino masses as M2 = 2000 GeV, 
Ml = 1000 GeV, {a) the slepton masses from the neutralino and chargino sector 
as Ml = 1500 GeV, and = 1000 GeV. 

As mentioned in the Introduction, we take into account of not very large vahies 
of tan/? (i.e. tan/3 < 10). For the purpose of showing the mass spectra, we shall 
take tan/3 = 2. To give feeling of the sensitivity of the mass on tan/3, we shall 
also consider tan/? = 10 regime. In the numerical analysis, we concentrate on two 
specific cases: in the first part, we carry out the analysis by letting A of a real 
parameter, then in the second part we take into account of the case for which A is 
complex. 
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3.1. The Case of Real A 

In the first part of our analysis, we take A as a real parameter {ipx = 0), and let 
all the other phases in the theory of being complex. Namely, fAt = 'PAb=fA^—Vi = 
(p2='-p. In Fig. 1, wc show the dependence of \eEDM\ on (p and A for selected values 
of the vacuum expectation value of the singlet (vg) when tan/3 = 2 (left panel) and 
tan/3 = 10 (right panel). In the figure, we obtain 3-dimensional surfaces for each 
value of Vs, which we choose within the v ^ Vs S interval. For instance, the top 
surface corresponds to Vg = 4u, whereas the bottom represents Vs = v. Here, ip 
changes from to 2n, and A from 0.12 to 0.82. As both panels of Fig. 1 suggest 
the upper bound on \eEDM\ gradually increases, with the increasing values of A, 
and of Vg. One notes that the present boimd of \eEDM\ is consistent with the 

present experimental bound in the full (p range, for all values v ^ Vs ^ 'iv, and 
0.12 < A < 0.82, at both tan/3 = 2 (left panel) and tan/3 = 10 (right panel) regimes, 
provided that the gaugino masses are of C?(TeV). 




Fig. 1. The dependence of \eEDM\ on ip and A for selected values of the vacuum expectation 
value of the singlet when tan (3 = 2 (left panel) and tan /3 = 10 (right panel) . 

The analysis of Fig. 1 gives a general idea of the parameter domain of the 
\eEDM\ and </?-A plane, when ipx = 0, and changes in its full range. With this 
input in mind, we choose two values of A (A = 0.12, and A = 0.45), and of Vs 
{vs ~ V and Vg = 3w), in the parameter space allowed by the eEDM constraint, 
corresponding to the low and high values of A and Vg, respectively, to analyze the 
dependence of the mass (m/i j at the tree and one- loop levels, and of the CP odd 
components of the lightest Higgs boson (h^) on ip at both tan 5 regimes. 

We would like to note that, in the remaining part of the analysis the variation 
of mh^ as well as ps and p5, with (f is displayed for several values of At taking into 
account of the strong dependence of the radiative corrections on the stop splitting, 
as will be seen in Eq. (22). 

Therefore, in each of the following plots, we first consider the case for which At 
and Vg are of comparable size, (i.e, At = Vg) when Vg = 175 GeV, and Vg = 525 GeV. 
Next, to determine how the increase in Ai affects thc^ radiative corrections, wc 
concentrate on two specific values of At corresponding to Ai = lOSOGeV. and 
At = UOOGeV for both Vg = 175 GeV and Vg = 525 GeV cases. 

In Fig. 2, we show the dependence of rrih^, (at tree and at one-loop levels) on p 
at tan (3 = 2 (left panel) and tan/3 = 10 (right panel), when Vg = v = 175 GeV, and 
A = 0.12, for selected values of At. In both panels of Fig. 2, the three upper curves 
with respect to the mid-point, from bottom to top, represent At = Vg = 175 GeV, 
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At = 1050 GcV, and At = 1400 GeV values of At respectively, whereas the lowest 

curve is for tlw^ trcc-lcvcl. 
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Fig. 2. The dependence of m/j^ on ip, for selected values of At, at tan /3 = 2 (left panel), and 
tan /3 = 10 (right panel) when Vs = 175 GeV and A = 0.12. Here, the three upper curves with 
respect to the mid-point, from bottom to top, are for At = Vg = 175 GeV, At = 1050 GeV, and 
At = 1400 GeV, respectively, whereas the lowest curve is for the tree-level. 



As Fig. 2 suggests when A = 0.12, and Vg = 175 GcV. to;,, grows from 
(42, 46, 54) to (69, 65, 56) GcV at tan /3 = 2, and from (71, 85, 109) to (135, 130, 115) 
GeV at tan/3 = 10, for the values of At = 1400 GeV (the top curve with respect 
to the mid-point), At = 1050 GeV (the second curve below the top curve), and 
At = Vs (the third curve below the top curve), respectively, as y ranges from to 
TT. On the other hand, remaining around 35 GeV for tan /3 = 2, grows up until 
87 GeV, at the tree level, for tan /3 = 10. It can also be observed ifrom Fig. 2 that 
TOftj is quite sensitive to the variations in </? for A = 0.12, and Vg = 175 GeV for 
both tan/3 = 2 (left panel) and tan/3 = 10 (right panel) regimes. One notes that 
the splittings between the tree and one-loop spectra are quite small. However, as 
Vg gets higher values, the strength of the radiative corrections are affected, which 
causes larger splittings between the tree and one-loop mass spectra. 
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Fig. 3. The dependence of rtifi^ on if, for selected values of At, at tan/3 = 2 (left panel), and 
tan/3 = 10 (right panel), when iis=525 GeV and A= 0.12. Here, the bottom, the middle and the 
top curves, are for At = Vs = 525 GeV, At = 1050 GeV, and At = 1400 GeV values, whereas the 
lowest curve is for the tree level. 



For instance, in Fig. 3, we show the dependence of to/i^ on ip, for Vg ~ iv = 
525 GeV, and A = 0.12. In Fig. 3, we select three specific values of At- namely 
At = 1400 GeV (the top curve). At = 1050 GeV (the second curve below the top 
curve). At = Vg GeV (the third curve below the top curve) Here, the lowest curve 
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represents m^^ at the tree-level. It can be seen from the left panel of Fig. 3 that as 
Vs increases the gap between the tree and the one-loop levels enlarges. For instance, 
being around at most ~ 20 GeV for Vg — v and At — Vg (left panel of Fig. 2), it 
rises to ~ 50 GeV for Vs = 3v at tan/3 = 2 (left panel of Fig. 3). This behaviour 
occurs also in the tan/3 = 10 regime for which case the gap between the tree and 
one-loop levels increases (right panel of Fig. 3). 

A comparative analysis of left panels of Fig. 2 and Fig. 3 shows that when At 
and Vg are of comparible size, the variation of ruh^ with (p is much more slower as 
compared to the cases for which At > Vg {At = 1050 GeV, and At = 1400 GeV). 
This saturation effect can be understood by observing that the radiative corrections 
depend strongly on the stop splitting A?, which depends explicitely on (pxt such that 



c5 = A?(7r)/A?(0) 



\At\-\X\vgCotp 1 



|A| cot 



\At\ + \X\Vg COt/3 



, |A|cot/3 



(23) 



This quantity particularly implies that the strength of the radiative corrections 
modify from one CP conserving point to the next. It can also be seen from (23) 
that when \ At \ and Vg, are of comparible size, as tan/3 increases, |A| cot /3 decreases 

and indeed, S approaches to unity in the high tan/3 limit. One notes that the tree 
level mass of the model is proportional to A^w^, whereas the CP violating entries 
of the radiative corrections of the Higgs mass-squared matrix grow with the term 
|j4i|usA. Naturally, with the increase in At, \At\vgX, term also increases. That is 
the radiative corrections enhance through the At term, which does not cause too 
big splittings among the radiative corrections, when Vg = v and A — 0.12. However, 
as Vg gets higher values, the seperation from the tree-level enlarges. 

In Fig. 4, we show the dependence of CP-odd components (ps and ps) of /15 
on (f, for A = 0.12, and Vg = v = 175 GeV, when tan/3 = 2 (left panel) and 
tan/3 = 10 (right panel). In both panels the upper and the lower curves represent 
the ps, and p5 components of ^,5, respectively. We select three values of At: At = 
Vg = 175 GeV, At = 1050 GeV, At = 1400 GeV, from top to bottom for the upper 
curves (pa), and from bottom to top for the lower curves (ps). 
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Fig. 4. The dependence of pa (the upper curves) and ps (the lower curves) on ip for Vs=l75 GeV, 
and A= 0.12, for selected values of At, when tan/3 = 2 (left panel), and tan /3 = 10 (right panel). 
Here, At = Vs = 175 GeV, At = 1050 GeV, At = 1400 GeV, from top to bottom for pa, and from 
bottom to top for ps. 

As the left panel of Fig. 4 suggests, when At = 1400 GeV, and A, = 1050 GeV 
/i5 has %98 and %96 components, respectively, at tan/3 = 2 and = 0. When 

= Vg, p3 decreases to %92. The left panel of the Fig. 4 also indicates that p5 
ranges in between %2 and %20, depending on the strength of Af On the other 
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hand, as wc observe from the right panel, when tan/3 = 10, rises nearly to the 
%100 line for At = 1050 GeV (and for At = 1400 GeV, as well), and correspondingly 
P5 remains in the vicinity of the %0 line at ip ^ 0. For higher values of the ps 
component of /15 increases, and this increase in the p5 is compensated by p3, as 
expected (see Fig. 5). 




0123456 0123456 

Fig. 5. The dependence of pa (the lower curves) and (the upper curves) on ip for Vs = 525 GeV, 
and A= 0.12, for selected values of At, when tan/3 = 2 (left panel) and tan/3 = 10 (right panel). 
Here, At = Vs = 525 GeV, At = 1050 GeV and At = 1400 GeV, from bottom to top for ps, and 
from top to bottom for p5, with respect to the mid-point. 



In Fig. 5, we show the dependence of p^, and ps components of /15 on ip, for 
A — 0.12, and Vg = 3v = 525 GeV, when tan (3 = 2 (left panel), and tan /3 = 10 (right 
panel). Here, the lower and the upper curves represent the p3, and ps components 
of respectively. Similar to the Fig. 4, in Fig. 5, we select three values of At: 
At = Vg, At = 1050 GeV, At = 1400 GeV, from bottom to top for the lower curves 
(ps), and from top to bottom for the upper curves (ps), with respect to the mid- 
point. It can be noted from the left panel of Fig. 5 that as Vs gets higher values the 
P5 component of /15 increases, with the decrease of its ps component, as compared 
to the previous case (see Fig. 4). For instance, when tan/3 = 2 (left panel of Fig. 
5), and At = 1400 GeV, ps decreases to ~ %63 whereas ps increases to ^ %37 at 
if = 0. A more spectacular side of Fig. 5 arises for tan /3 = 10 (the right panel) 
where one observes a very slow variation of <p as compared to the previous case (the 
right panel of Fig. 4). 

A comparative analysis of Figs. 4 and 5 suggest that for small A (A = 0.12) 
and higher Vg {vg = 3w), the ps component of increases, and this increase is 
compensated by the decrease in P3. Such kind of effect can be dominantly seen in 
the high tan/3 regime (see the right panel of Fig. 5). 

One notes that in the MSSM there is a single CP-odd component, and naturally 
the results of the CP-invariant theory are expected to be recovered at the CP- 
conserving points {(p = 0,7r, 27r) However, there are two CP odd components 
in the NMSSM (p3 and ps), and these two odd components mix each other at the 
CP conserving points, depending on the relative strengths of the other one-loop 
corrections. Therefore, one necessarily does not recover the results of the CP- 
invariant theory at the CP conserving points, due to the CP violating mixings of 
Ps and p5. 

Until now, we have carried out our analysis for small values of A, when (fx = 0, 
and all the other phases in the theory are assumed to be complex. In the following, 
we take into account of higher values of A. For this reason, we fix A = 0.45 and 
carry out the similar analysis for the second part of the work. 

In Fig. 6, we show the dependence of m^j_ (at tree and at one-loop levels) on ip for 
Vg = V = 175 GeV, and A = 0.45, when tan/3 = 2 (left panel) and tan/3 = 10 (right 
panel). We again select three values of At, namely At = 1400 GeV (the top curve 
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Fig. 6. The dependence of m/i^ on tp, for selected values of At, at tan/3 = 2 (left panel), and 
tan /3 = 10 (right panel) when Vs=^ and A = 0.45. Here, the three upper curves with respect to the 
mid-point, from bottom to top, are for At=Vs = 175 GeV, At = 1050 GeV, and At = 1400 GeV, 
respectively, whereas the lowest curve is for the tree level. 



with respect to the mid-point), At = 1050 GeV (the second curve below the top 
curve), At = Vs = 175 GeV (the third curve below the top curve). Here, the lowest 
curve is for the tree-level. It can be observed from Fig. 6 that m/i^ still remains 
sensitive to the variations in (f for higher values of A (A = 0.45) at tan/3 = 2. As 
tan/3 increases, the dependence of top squark masses on A, as well as Vg weakens, 
therefore these elements of the mass-squared matrix become more sensitive to the 
choice of these parameters. Hence, for tan/3 — 10 (right panel of Fig. 6), one 
observes a slow variation with (p, and quite large splitting between the tree and 
loop levels as compared to A = 0A2 case (see Fig. 2, right panel). Similar analysis 
can be performed for Vg ~ 2v, and A = 0.45, in which case one observes a weak 
dependence of (p on ruh^ at both tan/3 = 2 and tan/3 = 10 regimes (see Fig. 7). 
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Fig. 7. The dependence of m^^ on tp for Vs= 525 GcV, and A= 0.45, for selected values of At, 
when tan /3 = 2 (left panel), and tan /3 = 10 (right panel). Here, the bottom, middle and the top 
curves are for At = Vs = 525 GeV, At = 1050 GeV, and At = 1400 GeV, respectively, whereas 
the lowest curve is for the tree-level. 

One notes from the left panel Fig. 7 that the increase in Vg affects the radiative 

corrections, which causes quite large splittings between the tree and the loop levels, 
at tan/3 — 2. On the other hand, it is seen from the right panel that when tan/? = 
10, the variation of ip is much more flat, as compared to right panel of Fig. 6. 

In Fig. 8, we show the dependence of p3 and pf, on ip, for A = 0.45, and Vg = v, 
when tan /3 = 2 (left panel) and tan/3 — 10 (right panel). From the left panel it is 
seen that when At = 1050 GeV {At = 1400 GeV), the (f dependence of ps (the upper 
curves) and (the lower curves) is stronger as compared to At = Vg (the top curve 
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Fig. 8. The dependence of p3 and on (p for Vs=v, and A= 0.45, for selected values of At, when 
tan/3 = 2 (left panel), and taii/3 = 10 (right panel). The upper (lower) curves are for pz{pz) 
in the left panel, whereas they are for psips) in the right panel. Here, At = Vs = 175 GeV, 
At = 1050 GeV, At = 1400 GeV, from top to bottom for pa (ps), with respect to the mid-point. 



with respect to the mid-point). It is seen that as tan/3 increases, the contributions 
of the CP odd components change. For instance, at tan/3 — 2 (left panel), and 
At — Svg p3 starts from %98 at = and falls to %65 at (p = n (ps varies from %2 
to %37 in this interval). In passing to tan/3 = 10 regime (right panel) /15 is seen to 
gain non-negligible ps composition. For instance at At = 1400 GeV, ps starts with 
%50 at = 0, and rises to %90 at </? = tt. Similar to the previous observations, it 
can be observed that when At and Vg are of comparible size, the strengths of the 
CP violation mixing is weak and the variations of ps and ps with ip are quite slow. 
As the figure suggests, the higher At the larger the CP violating mixings between 
the CP odd components. One notes that similar analysis can be carried out for 
the higher values of Vs- It can be seen that, as Vs increases, ps component of hs 
increases, and this increase in the p5 component is compensated by pa. 

3.2. The Case of Complex A 

In the second part of the analysis, we take A as a complex parameter, with the 
phase ipx. Then, we let all the other soft phases in the theory of being zero {<PAt = 
<PAi,=<^Ae=¥'i= ¥'2=0). to explore the dependence of (p\ on various parameters. 




Fig. 9. The dependence of \eEDM\ on (fx and A for selected values of the vacuum expectation 
value of the singlet when tan j3 = 2 (left panel) and tan /3 = 10 (right panel) . 
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In Fig. 9, we show the dependence of \eEDM\ on (fx and A for selected vahies 
of Vs when tan/? = 2 (left panel) and tan/3 — 10 (right panel). Similar to the 
former analysis, we fix fc = 0.63, and consider all vahies of A in the 0.12 < A < 0.82, 
interval, as (p\ changes from to 27r. In both panels of the Figure, the top surface 
corresponds to Vg = Av, whereas the bottom represents Vg = v. A comparative look 
at both panels of Fig. 9 suggests that the present bounds of \eEDM\ is satisfied in 
all the parameter domain at both tan (3 regimes, when (p\ changes in its full range 
and A from 0.12 to 0.82, provided that the gaugino masses are of 0{TcV). 

In the following, we will perform a similar analysis that we have carried out for 
the real A case. To understand the effects of (p\ on m/i^, as well as and p^, wc first 
start with a low value of A, that is wc set A = 0.12, and again consider two specific 
values of Vs, namely Vs = 175 GeV, and Vg = 525 GeV. We show the dependence 
of ruh^ (at tree and at one- loop levels) on ipx, in Fig. 10 (when Vs = v), and Fig. 
11 {vg = 3v) at tan/3 = 2 (left panels) and tan/3 = 10 (right panels). 
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. 10. The dependence of mji^^ on tpx, for selected values of At, at tan,i3 = 2 (left panel) and 
f3 = 10 (right panel), when Vs = 175 GeV, v and A= 0.12. Here, the bottom, middle and 
top curves (with respect to the mid-point) are for At = Vs = 175 GeV, At = 1050 GeV, and 
= 1400 GeV, respectively, whereas the lowest curve is for the tree-level. 



As the left panel of Fig. 10 indicates m/i^ is quite sensitive to the variations 
of (fx when Vg = v, and tan/3 = 2. For tan/3 = 10 regime, it is seen that the 
variation of ipx with nih^ is slower at both tree and one-loop levels, particularly in 
the 1 < ~ 5 interval. 
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Fig. 11. The dependence of m/j^ on (px, for selected values of At, at tan/3 = 2 (left panel) and 
tan /? = 10 (right panel), when Vs = 525 GcV and A= 0.12. Here, the bottom, the middle and 
the top curves, are for At = = 525 GeV, At = 1050 GeV, and At = 1400 GeV, respectively, 
whereas the lowest curve is for the tree level. 



In passing to Vg = 3v case (see Fig. 11), although ruh^ remains quite flat for 
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tan/3 = 10 (right panel), one obtains a relatively fast variation for tan/3 — 2. It 
is seen that from the left panel of Fig. 11 that the peak value of at the one- 
loop level, as well as the tree level, is shifted through the maximal CP violation 
point at tan/3 = 2. A comparative analysis of Figs. 10 and 11 suggest that nih^ is 
quite sensitive to the variations in ipx for A = 0.12, when Vg = v and Vs — 3w at 
tan/3 = 2 (left panels) regime. Such kind of sensitivity weakens at tan/3 = 10 (right 
panels) for Vs ^ 2v. 

We show the dependence of CP-odd components (ps and ps), of /15 on (f\, for 
A = 0.12, and Vs = v in Fig. 12, and Vg = 3v in Fig. 13, respectively, when 
tan/3 = 2 (left panels) and tan/3 = 10 (right panels). As can be observed from the 
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Fig. 12. The dependence of pz (the upper curves) and (the lower curves) on ip\ for Vs = 
175 GeV, and A= 0.12, when tan (3 = 2 (left panel), and tan /3 = 10 (right panel). Here, At = Vs = 
175 GeV, At = 1050 GeV, At = 1400 GeV, from bottom to top for ps, and from top to bottom 
for p5, with respect to the mid-point. 



left panel of Fig. 12, when tan/3 = 2, component of on the average ranges 
from %98 to %90, as At changes from 1400 GeV to 175 GeV (from top to bottom 
for the upper lines) at ipx = 0. In passing to tan /3 = 10 regime, it is seen that /15 has 
nearly %100 ps component at <f\ = 0, as At changes from 1400 GeV to 1050 GeV 
(the first and the second curve from the top of the upper lines representing p.3). 
Moreover, is seen to gain non-negligible ps component in such a way that ps and 
P5 components are strongly mixed up in this regime (right panel). 
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Fig. 13. The dependence of pa (the lower curves) and pB (the upper curves) on (px for Vs = 
525 GeV, and X= 0.12, when tan/3 = 2 (left panel), and tan /3 = 10 (right panel). Here, At = Vs = 
525 GeV, At = 1050 GeV, At = 1400 GeV, from bottom to top for ps, and from top to bottom 
for p5, with respect to the mid-point. 



On the other, considering the higher values of Vg (see Fig. 13), one notes that 
the p5 component of increases, and this increase in the p5 is compensated by 
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P3 when tan/3 = 2 (left panel). For tan/3 = 10 (right panel) regime, it can be 
observed that p^, component of has a slow variation in the vicinity of %100 line, 
whereas pz remains around %0 line, for each value of At, changing from 525 GeV 
to At = 1400 GeV. 

In the last part of our analysis, we fix A = 0.45, and carry out the similar analysis 
for Vs = V and Vs = 3u. We show the dependence of ruh^ (at tree and at one- loop 
levels) on ^px for Vg = v m Fig. 14, and Vs = 3t) in Fig. 15, when tan/3 = 2 (left 
panels) and tan/3 = 10 (right panels). In Figs. 14 and 15, At = 1400 GeV (the top 
curve with respect to the mid-point), At = 1050 GeV (the second curve below the 
top curve), whereas the third curve below the top curve is for At = = 175 GeV 
in Fig. 14, and At = Vs = 525 GeV in Fig. 15. In both Figures, the lowest curve 
represents m/,^ at the tree-level. 
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Fig. 14. The dependence of mji^ on ipx for Vs= v, and A= 0.45, for selected values of At, 
when tan /3 = 2 (left panel) , and tan ,8 = 10 (right panel) . Here, the bottom, the middle and 
the top curves, with respect to the mid-point, are for At = Vs = 175 GeV, At = 1050 GeV, and 
At = 1400 GeV values, whereas the lowest curve is for the tree level. 



200 I 1 — 5 — I 1 

180 - ''-J 

160 - / *'\ 

^ 140 - / / 

2. 120 - :•' / \ . ,.*' / 

E 100 - / \ " /' 

80 - / \ / 

60 :v' 



0123456 0123456 

Fig. 15. The dependence of mf^^ on (fix for Vs=3 v, and A= 0.45, for selected values of At, 
when tan/3 = 2 (left panel), and tan /3 = 10 (right panel). Here, the bottom, the middle and 
the top curves, with respect to the mid-point, are for At = Vs = 525 GeV, At = 1050 GeV, and 
At = 1400 GeV values, whereas the lowest curve is for the tree level. 

One notes from the left panels of Figs. 14 and 15 that as A increases (A = 0.45), 

one obtains a quite fast variation of with ipx both at the tree and loop levels, 
particularly in the low-tan /3 regime. The increase in A, as well as Vs affects the 
spectrum of /is in the sense that the shifting of the peak value of m/i^ (at the tree 
and the loop levels) towards the; maximal CP violation point is sharper as compared 
to the A — 0.12 case (see Fig. 11). Although the sensitivity much more weakens in 
the high tan /3 regime, such kind of shifting effect of can also be observed from 
the right panel of Fig. 15. 
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Finally, in Fig. 16, wc show the dependence of CP-odd components (p3 and p^) 
of on (fx, for A = 0.45, and Vg — v at tan/3 — 2 (left panel) and tan/3 = 10 (right 
panel). As can be observed from Fig. 16, the increase in A affects the CP odd 
components in such a way that ps and p5 mix significantly in the tan (3 = 2 regime. 
In passing to tan/3 = 10 regime, it is seen that such mixings weaken, however /15 is 
seen to gain non-negligible ps component as At increases. On the other hand, as Vg 
increases, it can be seen that p3 and p5 have similar ip\ dependencies for all values 
of At, which essentially follows from the dominance of the AfXvg term. 
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Fig. 16. The dependence of pa and ps on ipx for Vs= v, and A= 0.45, for selected values of At, 
when tan /3 = 2 (left panel), and tan/3 = 10 (right panel). The upper (lower) curves arc for P3(p5) 
in the left panel, whereas they are for p5( pa) in the right panel. Here, At = Vs = 175 GeV, 
At = 1050 GeV, At = 1400 GeV, from top to bottom for pa (ps), with respect to the mid-point. 



Before concluding, we would like to close this section with a brief discussion of 

the measurement of CP violating effects at colliders: Recently, the authors of Ref. 
[44] have presented a general formalism, the resonant CP violation, for analyzing the 
CP violating phenomena in the MSSM at high-energy colliders. Their formalism, 
which is developed from Ref. [45], can bo applied to models with an extended 
CP violating Higgs sector, including the MSSM with radiative Higgs-scctor CP 
violation. Resonant CP violation has not been studied in the NMSSM yet. The 
measurement of resonant CP violation effects can provide a useful tool for the 
determination of the CP odd compositions of the Higgs bosons in the NMSSM. 



4. Conclusion 

We study the explicit CP violation of the Higgs sector in the next-to-minimal 
supersymmetric model with a gauge singlet Higgs field. Our general discussion 
followed by the numerical estimates for various parameter planes show that: 

(i)The eEDM lies around the present experimental upper limits, provided that 
the gaugino masses are of ©(TeV). 

(M)When A is real, and all the other phases in the theory are complex, nih^ is 
quite sensitive to particularly at small A and Vs (A = 0.12, and Vg = v). However 
as Vs and tan/3 increases (for instance, A = 0.12, Vg = 3u and tan/3 = 10), the 
radiative corrections which are sensitive to variations in ip are suppressed. Indeed, 
m/j5 nearly remains constant for tan/3 = 10 at Vg = "iv. Considering the ps and p5 
components of /15 at small A and Vs (A = 0.12, and Vs = v), as At increases, the 
p:j component of h-, decreases, whereas its ps component increases. However, this 
increase in p5 can not be larger than %20, at tan /3 = 2, even for At = 1400 GeV. 
In passing to tan /3 = 10 regime, it is seen that the p3 component of /i5 tend to 
decrease more rapidly, whereas this decrease is compensated by the increase in ps, 
in the sense that the two components mix aromid = ti, when A = 0.12. On the 
other hand, for higher values of Vg [vg = 3ti), it is seen that the variation of the CP 
odd components is quite slow, as compared to the former case. 
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(in) When A is complex, and all the other CP violating parameters in the theory 
are chosen to be real, rrihr^ grows significantly as compared to the former case. It is 
observed that (px strongly affect the the tree and one-loop Higgs masses. Clearly, 
ipx is not only a physical phase which forms the source of CP violation at the tree 
level, but it also affects the one- loop radiative corrections via (p\t. Therefore, even 
if all the other phases are set equal to zero, the effects of ipx can be seen on the mass 
and the CP components of /15. The CP violating effects are particularly enhanced, 
as Vs increases. 

(w)Being the most economic extension of the MSSM, the next-to-minimal su- 

persymmetric model not only leads us to a wealth of CP violation opportunities, 
but it also offers a rich phenomenology for future colliders. 
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